EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-PH-EP/2004-0 1 3 
5th April 2004 

Measurement of the Branching Ratio and Form Factors for 
the Decay K L —> 7t ± 7t°e =F 'Ve('Ve) 



§ ' NA48 Collaboration 

X 



X 



J.R. Batley, R.S. Dosanjh, T.J. Gershon, G.E. Kalmus, C. Lazzeroni, D.J. Munday, E. Olaiya, 
M. Patel, M.A. Parker, TO. White, SA. Wotton 
Cavendish Laboratory, University of Cambridge, Cambridge, CB3 OHE, U.K. 1 

R. Arcidiacono, G. Bocquet, A. Ceccucci, T. Cuhadar-Donszelmann, D. Cundy 2 , N. Doble 3 , 
V. Falaleev, L. Gatignon, A. Gonidec, B. Gorini, P. Grafstrom, W. Kubischta, I. Mikulec 4 , 
q ! A. Norton, S. Palestini, H. Wahl 5 

ID i CERN, CH-121 1 Geneve 23, Switzerland. 

C. Cheshkov", P. Hristov", V. Kekelidze, L. Litov", D. Madigozhin, N. Molokanova, 

Yu. Potrebenikov, A. Zinchenko 
Joint Institute for Nuclear Research, Dubna, 141980, Russian Federation 

j^; P. Rubin 7 , R. Sacco 8 , A. Walker 

Department of Physics and Astronomy, University of Edinburgh, Edinburgh, EH9 3JZ, U.K. 1 

D. Bettoni, R. Calabrese, P. Dalpiaz, J. Duclos, PL. Frabetti 9 , A. Gianoli, M. Martini, 
L. Masetti 10 , F. Petrucci, M. Savrie, M. Scarpa 

Dipartimento di Fisica dell'Universita e Sezione delVINFN di Ferrara, 1-44100 Ferrara, Italy 

A. Bizzeti 11 , M. Calvetti, G. Collazuol 3 , E. Iacopini, M. Lenti, F. Martelli 12 , G. Ruggiero, 

M. Veltri 12 

Dipartimento di Fisica dell'Universita e Sezione delVINFN di Firenze, 1-50125 Firenze, Italy 

K. Eppard, M. Eppard 6 , A. Hirstius 6 , K. Kleinknecht, U. Koch, L. Kopke, P. Lopes da Silva, 
P. Marouelli, I. Mestvirishvili, C. Morales, I. Pellmann 13 , A. Peters 6 , B. Renk, S.A. Schmidt, 

V. Schonharting, R. Wanke, A. Winhart 
Institut fur Physik, Universitat Mainz, D-55099 Mainz, Germany 15 

J.C. Chollet, L. Fayard, G. Graziani, L. Iconomidou-Fayard, G. Unal, I. Wingerter-Seez 
Laboratoire de V Accelerateur Lineaire, IN2P3-CNRS, Universite de Paris-Sud, F-91898 

Orsay, France 16 

G. Anzivino, P. Cenci, E. Imbergamo, G. Lamanna, P. Lubrano, A. Mestvirishvili, A. Michetti, 

A. Nappi, M. Pepe, M. Piccini, M. Valdata-Nappi 
Dipartimento di Fisica dell'Universita e Sezione delVINFN di Perugia, 1-06100 Perugia, Italy 

R. Casali, C. Cerri, M. Cirilli 6 , F. Costantini, R. Fantechi, L. Fiorini, S. Giudici, I. Mannelli, 

G. Pierazzini, M. Sozzi 



Dipartimento di Fisica, Scuola Normale Superiors e Sezione dell'INFN di Pisa, 1-56100 Pisa, 

Italy 

J.B. Cheze, M. De Beer, P. Debu, F. Derue, A. Formica, G. Gouge, G. Marel, E. Mazzucato, 

B. Peyaud, R. Turlay 14 , B. Vallage 
DSM/DAPNIA - CEA Saclay, ¥-91191 Gif-sur-Yvette, France 

M. Holder, A. Maier, M. Zioikowski 
Fachbereich Physik, Universitat Siegen, D-57068 Siegen, Germany 17 

C. Biino, N. Cartiglia, M. Clemencic, F Marchetto, E. Menichetti, N. Pastrone 
Dipartimento di Fisica Sperimentale dell'Universita e Sezione delVINFN di Torino, 1-10125 

Torino, Italy 

J. Nassalski, E. Rondio, W. Wislicki, S. Wronka 
Soltan Institute for Nuclear Studies, Laboratory for High Energy Physics, PL-00-681 Warsaw, 

Poland 1 * 

H. Dibon, M. Jeitler, M. Markytan, G. Neuhofer, M. Pernicka, A. Taurok, L. Widhalm 
Osterreichische Akademie der Wissenschaften, Institut fur Hochenergiephysik, A- 1050 Wien, 

Austria 19 



Submitted to Physics Letters B 



Abstract 



The K L -> 7t ± 7T°e T 'Ve( , Ve) decay was investigated with the NA48 detector at CERN SPS 
using abeam of long-lived neutral kaons. The branching ratio Bt(Kl — > 7T ± 7t°e =F "v e ("v e )) = 
(5.21 ± 0.07 sta t ± 0.09 syst ) x 1 0~ 5 was fixed from a sample of 5464 events with 62 back- 
ground events. The form factors f s , f p , A g and h were found to be in agreement with 
previous measurements but with higher accuracy. The coupling parameter of the chiral 
Lagrangian L3 = (—4.1 ± 0.2) x 10~ 3 was evaluated from the data. 



1 Introduction 

The decay Kl — > nncv, called K e 4, is recognized as a good test for chiral perturba- 
tion theory (CHPT) and its predictions for long-distance meson interactions. In particular, it 
is used to determine the nn partial wave expansion parameters: threshold parameters, slopes 
and scattering lengths, where the S-wave nn scattering lengths can be further related to the 
quark condensate Q- The complete set of CHPT parameters has been calculated in the one- 
loop approximation C(p 4 ) and the form factors F and G and quark condensates in the two-loop 
approximation O (p 6 ) Q. 

Following the initial observation of charged K e 4 0, the process K + — > 7r + 7r~e + ~v e , called 
Kg 4 , was measured in ref. JU based on an event sample of 30,000 events and, more recently, a 
high-statistics experiment [5] detected 400, 000 such decays. Those experiments determined the 
K+ 4 decay rate, four form factors and the difference of the s- and p-wave phase shifts 5^ — 5] as a 
function of the mass of the pion pair M^. By fitting the Roy model (6| to the M^-dependence 
of phase shifts and assuming time-reversal invariance, they also evaluated the scattering length 
ag. 

After a low-statistics observation of the neutral K e4 decay K L — > 7r ± 7r°e =F 'V e (~v e ) Q, a 
more complete analysis was performed in ref. ID, where a sample of 729 events was used to 
determine the branching ratio, the threshold value g(M 7T7T = 0) of the g form factor, the relative 
form factors f s = f s /g, f p = f p /g and h. = h/g, and the M^-dependence of g. The neutral 
K e4 decay is well-suited for measuring the G form factor and the L 3 parameter of CHPT. 

This paper reports on the measurement of both the branching ratio and the form factors 
of neutral K e4 decays by the NA48 Collaboration at CERN, using a significantly larger data 
sample than previous measurements. In addition, the coefficient L3 of the chiral Lagrangian, 
sensitive to the gluon condensate, is evaluated with high accuracy. 
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2 Kinematics and parametrization of the decay cross section 

The matrix element for the decay is assumed to factorize into a leptonic term, describing 
the coupling of the W boson to leptons, and an hadronic term, accounting for hadronization of 
quarks into pions and representing the V — A structure iTTOl 

M = sin e c (7r7r|A A + V a |K)iWy a (1 - y 5 K, (1) 

where Gp is the Fermi weak coupling constant and C the Cabibbo angle. The vector compo- 
nent (7T7r|V|K) is parametrized in terms of one form factor H 

(7m|V A |K> = -LH £ A ^(p K y P7tl + p 7t2 )„(p 7tl - V n 2 ) P , (2) 
m K 

and the axial-vector part (7ra|A|K) in terms of three form factors: F, G and R 

(7t7t|A A |K) = [F( P7tl + Vn2 ) + G(p ni -p^J+R^K-Pn, -Pn 2 )] A , (3) 

m K 

where the R term is suppressed by the squared ratio of the electron mass to the kaon mass and 
can therefore be neglected. 

The differential cross section for the K e 4 decay was proposed El to be analysed in 
terms of five Cabibbo-Maksymowicz (C-M) variables: invariant masses of the dipion M.^ and 
the dilepton M. ev , the polar angles 0^ (9 e ) between the charged pion (electron) momentum in 
the nn (ev) centre-of-momentum frame and the dipion (dilepton) momentum in the kaon rest 
frame, and the azimuthal angle <b from the nn plane to the ev plane (Fig.[T]). 




Figure 1: Definition of the Cabibbo-Maksymowicz kinematic variables [9| used for the analysis 
of K e 4 decays. The angles Q n (9 e ) are between the charged pion (electron) momentum in the 
dipion (dilepton) centre-of-momentum frame and the dipion (dilepton) momentum in the kaon 
rest frame. The directed angle (b is from the nn plane to the ev plane. 



The 9^-dependence of the form-factors is made explicit by using a partial-wave expansion 
of the hadronic matrix element with respect to the angular momentum of the pion pair and 
restricting this expansion to s and p waves due to the limited phase space available in the K e 4: 

F = f s e i6s + f p e i5p cos Q n 
G = ge i6p 

H = he i& ? (4) 



2 



The G and H expansions contain only the p wave due to their antisymmetry with respect to 
pion exchange. Using the partial wave decomposition ©, an explicit expansion of the decay 
cross section in terms of the form factors f s , f p , g, H and 5 = 6 S — 6 P , and C-M variables M^, 
Mev, cos Q n , cos e and cb, was taken from ref. 0. The possible M^-dependence of g was 
accounted for by parameterizing g(M 7r7T ) = g(0)[l + A g (M^ 7t /4m^ : — 1 )], where A g , together 
with the form factors, has to be determined from a fit to the data, and m n stands for the average 
of the charged and neutral pion masses. 

3 The beam 

The NA48 experiment used for this investigation a 400 GeV/c proton beam from the 
CERN Super Proton Synchroton with a nominal intensity of 1 .5 x 1 12 protons per spill, deliv- 
ered every 16.8 s in 4.8 s long spills 1111 . Two kaon beams, one providing Kl decays and called 
the Kl beam, and another one, providing Ks decays, and called the Ks beam, were produced si- 
multaneously on two separate targets. For the K e 4 measurement only the Kl beam was relevant. 
The Kl beryllium target was located 126 m before the decay region. Charged particles were 
swept by dipole magnets, and the remaining neutral beam was defined by a set of collimators. 
The total flux of Kl's at the entrance of the fiducial decay volume was 2 x 1 7 per spill. 

4 The detector 

The detector system, located 114 m after the Ks target and extending 35 m downstream, 
consisted of two principal subsystems: a magnetic spectrometer and a spectrometer for neu- 
tral decays. In addition, there were scintillating hodoscopes, a hadron calorimeter, muon veto 
counters, beam veto counters, and a tagging station on the Ks beamline. 

The magnetic spectrometer was contained in a helium tank and consisted of a dipole 
magnet with a transverse momentum kick of 265 MeV/c and four drift chambers, each equipped 
with eight sensitive planes, arranged two before and two after the magnet. The momentum 
resolution of this spectrometer was between 0.5 % and 1 %, depending on the momentum, and 
the average plane efficiency exceeded 99 %. 

A scintillating hodoscope, consisting of two orthogonal planes of scintillating strips (hor- 
izontal and vertical) had a time resolution of 150 ps. Signals from quadrants were logically 
combined and used for triggering charged events in the first level trigger. 

An iron-scintillator hadron calorimeter, 6.7 nuclear interactions thick and located down- 
stream of both spectrometers, provided a total energy measurement, complementary to the elec- 
tromagnetic calorimeter. 

Muon veto counters, situated behind the hadron calorimeter, provided time information 
used to identify muons and to suppress backgrounds both in low-level triggers and offline. 

The fiducial decay region was surrounded by seven sets of iron-plastic veto scintillators, 
called AKL, used for identification of photons escaping this volume. 

The spectrometer for neutral decays consisted of a quasi-homogeneous ionization cham- 
ber calorimeter filled with 10 m 3 of liquid Krypton. Its length, amounting to 27 radiation lengths 
with a Moliere radius of 4.7 cm, ensured full containment of electromagnetic shower of energies 
up to 100 GeV, excluding detector regions close to the edges. The calorimeter was divided into 
13,212 cells, 2 cm x 2 cm transversally to the beam, read out individually. This calorimeter pro- 
vided the resolution of reconstructed energy cr ( E)/E = 9%/E©3.2%/v / E©0.42% and good 
reconstruction of the neutral vertex position along the beam. Signals from the calorimeter were 
digitized asynchronously by a 40 MHz flash ADCs and read out with online zero-suppression. 

A more detailed description of the apparatus can be found in ref. lfT2l . 
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5 The trigger 

Data were taken using the minimum-bias trigger ETOT, requiring a minimal energy de- 
posit of 35 GeV in the calorimeters, hit multiplicity in the first drift chamber, and a coincidence 
between opposite quadrants of the scintillator hodoscope. Since this trigger was downscaled by 
a factor of 30, a dedicated trigger KE4 was added to enhance statistics. The KE4 trigger used 
the neutral trigger system of NA48, which gave information about x and y projections of the 
energy deposited in the electromagnetic calorimeter. The requirement was at least 3, and not 
more than 5, clusters in any of the two projections, corresponding to the two photons from a n° 
decay and two charged particles, allowing for the loss of one cluster due to overlapping. This 
trigger was downscaled by a factor of 50. In the liquid crypton calorimeter readout the KE4 
trigger induced lower threshold than the ETOT trigger. In case of both triggers conditions were 
fulfilled, the lowest one was used. Therefore a small number of events giving both triggers were 
included to the KE4 sample. 

The minimum-bias trigger ETOT was assumed to be fully efficient. The efficiency of the 
KE4 trigger was measured relative to the ETOT trigger and found to be [98.76 ± 0.21 ) % for 
K e 4 events and (98.92 ± 0.01 ) % for events. The latter were used for normalization. 

6 Data sample 

The analysis described in this paper refers to the data collected in 2001. The sample of 
K e 4 events was selected, from both the KE4 and ETOT triggers, by applying the following cuts 
to reconstructed events which fulfilled the triggering conditions: 

1. Two well-reconstructed tracks of opposite charges. 

2. Four reconstructed clusters in the liquid Krypton electromagnetic calorimeter. 

3. Two photon clusters, each of energy between 3 and 100 GeV, not associated with the 
charged tracks. 

4. Tracks impacting the Krypton calorimeter between 15 cm and 120 cm from the beam 
axis. This cut eliminates clusters close to the edge and thus not fully contained in the 
calorimeter. 

5. A minimum distance between photon and the charged pion clusters of 15 cm, thus ensuring 
cleanliness of cluster reconstruction. 

6. A minimum distance of 5 cm between a photon cluster and the extrapolation of the charged 
tracks from before the magnet, partially removing background from decays K L — > 7t ± e T ^ e 
with two additional photons, at least one of them coming either from internal or external 
bremsstrahlung (K e 3+2Y)- 

7. The total energy deposited in the Krypton calorimeter had to be larger than 30 GeV. 

8. The energy over momentum ratio (E/p) for the electron candidate had to be larger than 
0.9 and smaller than 1.1 and for the charged pion candidate smaller than 0.8 (cf. fig. EI). 

9. A Xsn variable for the K L — > n + n~n° hypothesis was defined as 



with the invariant mass under the 3n hypothesis and the transverse momentum pj (cf. 
fig. Ell. Mk is the kaon mass, pj = 0.006 GeV/c is the modal value of the pj distribution, 
cr M = 0.0025 GeV/c 2 and <y p = 0.007 GeV/c. The cut xl n > 1 6 suppresses most of the 
background, where one of the charged pions is misidentified as the electron. 
10. The invariant mass of the two-photon system, at the vertex defined by the two charged 
tracks, had to be between 0. 1 1 and 0. 15 GeV/c, which ensures that the photons come from 
a 7i° decay. 




(5) 
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Figure 2: Distributions of E/p for data. The wide distribution on the left corresponds to pions 
from cleanly selected sample, and the narrow one on the right to electrons from the K e3 
sample. These data were used to train the neural network. The cuts on E/p are also shown. 



1 1 . The ratio Pt/E-v had to be between and 0.02, where p T is the total transverse momentum 
of all visible particles (two pions and an electron) and E v is the energy of the neutrino in 
the laboratory frame. This cut suppresses the K e 3 + 2 y background with one or both photons 
coming from accidental coincidence. In this case the energy taken by the photons may be 
large enough to lead to a negative E v (cf. fig. H|). The K e3+2Y event sample in fig. |U was 
selected using a neural network algorithm. The cut is not efficient for K e3+ 2 Y background 
events with only bremsstrahlung photons, which are rejected by cuts 6 and 10. 

Cuts were also made on the maximum time difference between calorimeter clusters be- 
longing to the same event and between clusters and tracks. Cluster quality criteria were met 
as were requirements for spatial cluster separation, cluster versus track spatial matching, and 
vertex position and quality. Cuts 6, 10 and 11 above suppress background from K e 3 +2Y down 
to the level of 1.5% and 2.2% for the KE4 and ETOT triggers, respectively. This was estimated 
from Monte Carlo by normalizing the Pt/E v spectrum from the K e3+ 2 T to the one from K e4 in 
the region of Pt/E^, < and calculating the contamination for Pt/E-v > 0. 

Requirements 8 and 9 eliminate most of the background from the channel. The effects 
of the xin an d E/P cuts are illustrated in figs[5]where the background is shown as a function 
of E/p with xln > 1 6 (left) and as a function of 1 lx\ n (right). 

In order to diminish this background further, a neural network algorithm was applied 
[ 1 3 1 . A 3-layer neural network was trained on cleanly selected K e 3 and data samples to 
distinguish pion and electron electromagnetic showers in the liquid Krypton calorimeter. Both 
the and the K e 3 samples were taken during the same run period as the signal events. The 
algorithm used geometric characteristics of showers and tracks and E/p of tracks on the input 
and returned a control variable which was around for pions and around 1 for electrons. 

The background from was estimated as a function of x\n by extrapolating the pion 
tail shape under the electron E/p peak as determined from data. This was done without the 
use of the neural network, since the background would be too small otherwise. Then, with the 
neural network operating, the previously determined dependence on xl n wa s fit to the data. The 
background of charged pions misidentified as electrons in a pure sample amounted to 1.2% 
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Figure 3: Distribution of measured mass of three visible particles, assuming the 3n hypothesis, 
versus their total p T , for Monte Carlo K e4 and events. The ellipse defines the cut xj, n > 1 6 
which distinguishes K e4 events (outside the ellipse) from the background located around 
the kaon mass (inside the ellipse). 
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Figure 4: Distributions of p-p/tv f° r the accepted and reconstructed Monte Carlo K e4 events 
(solid) and the K e3+ 2y background events from data (dashed). For K e4 , resolution smearing may 
occasionally lead to events with small negative values for E^, and large negative pj/E^,, which 
are hardly seen in the plot. The sharp edge of both distributions at pj/E^, = is due to the 
kinematic suppression of large values of E^,. 
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Figure 5: The background shown as functions of the E/p with xin > ^ (l e ft) an d of the 
1 /xItt cut (right). For the E/p distributions (left), the data were corrected for the downscaling 
factor when E/p < 0.7 which results in larger errors in that domain. In the right-hand figure, 
the lower curves show the remaining K e 3 and backgrounds, the middle curve corresponds 
to the remaining signal, and the upper one shows all events. Both plots show data. 



and 1.1% for the KE4 and ETOT triggers, respectively. 

After applying all selection criteria, the data samples amounted to 2089 and 3375 events 
for the KE4 and ETOT triggers, respectively. The estimated numbers of background events were 
25 ± 10 and 37 ± 5, respectively. In total, the sample of K e 4 events amounted to 5464 events 
with 62 background events. 



7 Monte Carlo simulation 

Monte Carlo simulations were used for: 

- Determination of the acceptance for evaluation of the K e 4 branching ratio and form factors, 

- Estimation of the background, 

- Estimation of radiative corrections. 

The simulation code is based on event generators for neutral kaon decay channels and full 
GEANT simulation ifLU of all electromagnetic processes in the NA48 detector, including cas- 
cades in the calorimeter. Samples generated included: more than 1 million K e 4 events using 
the form of the decay matrix element as calculated by Pais and Treiman ifTOl. 0.5 million 
K e 3+2y events where one photon was from inner bremsstrahlung and another one from the 
bremsstrahlung of a charged particle in the detector's material, and 0.5 million events. 
Wherever possible, the data were used for background studies in order to be independent of 
Monte Carlo simulation. 



8 Branching ratio 

The branching ratio of the K e4 channel was determined by normalizing it to the 
channel: 

„ N(K e4 ) Q(K^) £(K^) 

Br(K -> = ^ ■ TiK^T ^ (6) 
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where N stands for the overall number of accepted events, properly corrected for downscaling, 
a for acceptance, and e for trigger efficiency. For the branching ratio of the reference channel, 
the value Br(K 7t3 ) = (12.58 ±0.1 9) % was used lfT5l . The reference sample of was selected 
using similar cuts 1-11 as for the K e4 but requiring xl n < 5 in cut 9 and no cut on E/p. 

The acceptances, as calculated from the ratios of accepted to generated Monte Carlo 
events, are equal to (3.610 ± 0.017) % and (5.552 ± 0.033) % for the K e4 and channels, 
respectively. 

From this, the branching ratios for the KE4 and ETOT trigger samples were found to be 
(5.30 ±0.12 stQt ± 0.11 syst) x 10~ 5 and (5.15 ± 0.09 stat ± 0.1 2 syst ) x 1 0~ 5 , respectively, and 
the overall branching ratio is equal to 

Br(K e4 ) = (5.21 ± 0.07 stat ± 0.09 syst ) x lO" 5 (7) 

where contributions to the systematic error, in units of 10~ 5 , are as follows: 

branching ratio: 0.079 

K e4 form factors: 0.021 

Background from K e 3 +2y : 0.019 

Background from K e4+y : 0.01 1 

Monte Carlo statistics for K e4 : 0.024 

Monte Carlo statistics for K^: 0.030 

Trigger efficiencies: 0.005 

Background from K^: 0.001 

These systematic errors for branching ratios for the KE4 and ETOT triggers were deter- 
mined with independent sets of data and Monte Carlo. The branching ratio includes radiative 
events K L — > n°n ± e T ^ e [y e )y (called K e4+Y ) left in the sample after all cuts. Their contribu- 
tion is accounted for in the systematic error and was estimated using the ratio of accepted K e4 
to K e 4+y, known from the Monte Carlo radiative event generator PHOTOS 11160 . The ratio of 
decay rates can be calculated using formulae of ref. ffol. It was found that the fraction of all 
K e4+y events in the final sample, including those with the radiative photon undetected due to 
acceptance or cuts, is 0.89 ± 0.02 %, in agreement with numbers cited in ref. Q. For the overall 
sample this corresponds to 48 events or a 0.042 x 1 0~ 5 contribution to the systematic error. The 
systematic uncertainty due to radiative corrections was estimated as ±25% of the maximal con- 
tribution to the branching ratio coming from the remaining K e4+Y background. This uncertainty 
represents an upper bound based on our experience with other decays. 

Clearly, the systematic uncertainty is dominated by the error in the branching ratio. 

Consistent results were obtained in two, independent analyses. 

9 The form factors 

The form factors were estimated by fitting the differential distributions, as found in the 
analytic expression in ref. [4], to the empirical distributions of the C-M variables. The detector 
acceptance was accounted for by multiplying the theoretical function by the distribution of 
accepted events generated by Monte Carlo with flat form factors. In order to account for K e4+Y 
events, the radiative generator PHOTOS ifToll was used. We generated radiative events with one 
bremsstrahlung photon emitted by a charged particle in the final state. 

Since the data sample of K e4 events was not large enough to allow a five-dimensional 
analysis, a simultaneous fit to all one-dimensional projections was performed. It was found, and 
checked with Monte Carlo, that the maximum correlation coefficient between points on projec- 
tions depends on the number of bins n as 0(1 /n 2 ) and therefore can be neglected for n > 1 0. 
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Figure 6: Distributions of the Cabibbo-Maksymowicz variables: (upper left), M CT (upper 
right), cos 9^ (middle left), cos e (middle right) and <b (bottom) for data from both KE4 and 
ETOT triggers (points with error bars) with fits (histograms). Acceptance was accounted for in 
the fit. 
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The theoretical curves were found for each projection by integrating over the remaining four 
variables. The method was checked by fitting the two-dimensional distributions of M. nn vs. 
M.^ and Q n vs. 9 e , after integration over the remaining three variables, and finding results in 
satisfactory agreement with a simultaneous one-dimensional fit. The one-dimensional distribu- 
tions of the C-M variables, experimental and fits, are presented in figs|6j The values of the form 
factors were found from the fit to data from combined KE4 and ETOT triggers to be: 

f s = 0.052 ± 0.006 stat ± 0.002 syst 
f p = -0.051 ±0.011 stat ±0.005 syst 
A g = 0.087 ±0.0]9 stat ±0.006 syst 

h = -0.32±0.12 stQt ±0.07 syst (8) 

withx 2 /ndf = 137/146. 

The phases 6 S and 5 P were not determined in this analysis. The factor g(0) can be deter- 
mined in a model-dependent way, using the branching ratio and predictions of CHPT lfT9ll . 

The systematic errors are dominated by the background, with minor contributions from 
Monte Carlo statistics and the shapes of the E/p and p T distributions used for background sub- 
tractions. As with the branching ratio, the contribution from the radiative K e4+Y events to the 
systematic errors was estimated as ±25% of the difference between the form factors calculated 
with and without radiative events in the Monte Carlo background. As a cross check, we fitted 
separately the KE4 and ETOT trigger samples obtaining results consistent within the statisti- 
cal errors. Additional checks were performed by inspecting ratios of C-M variable distributions 
between trigger samples, and between data and MC weighted by fitted form factors, and no sig- 
nificant discrepancies were found. The stability of the results was also examined by increasing 
and decreasing the amount of the background through varying the Xs n cut an( ^ f° un d to 
be satisfactory. Also the fitting procedure was checked against background addition. Radiative 
corrections affect the values of form factors within one standard deviation of the statistical error. 

Assuming the K e4 hypothesis, there are two kinematically allowed solutions for the kaon 
energy and hence the neutrino energy. The choice has no effect on but can slightly affect 
other C-M variables. Both solutions with equal weights of 0.5 were used for the form factor fits. 
The systematic effect of making the wrong choice was examined using Monte Carlo and found 
to be negligible. 

All steps of the analysis, viz. event selections, backgrounds and fits, were done twice and 
independently, and the results were in good agreement. 

10 Discussion and conclusions 

The K e4 branching ratio measured by NA48 is consistent with previous measurements 
000 within errors and is more accurate by a factor of 2.5 than that of ref. 0> both statistically 
and systematically (cf. fig. 0). The form factors f p , A g and h, have also significantly higher 
accuracy and agree within errors with ref. ID, whereas the value of f s differs by two standard 
deviations. 

We found a non-zero value of f s , allowing for the violation of the AI = 1 /2 rule at the 
percent level. As discussed in ref. admixtures of AI = 3/2 and AI = 5/2 to the K e4 decay 
amplitude, at the level comparable to that of the K — > nn decays, can be expected. 

We note good agreement of our h. value with previous neutral and charged K e4 studies 
BJI51IE1 and with the theoretical prediction 0, essentially independent of the coefficients of 
the chiral Lagrangian. 
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Figure 7: K e4 branching ratios from refs J7J and and from the present work (results from 
data using both triggers are shown). The shaded belt corresponds to the overall NA48 result. 



The K e4 decay is helpful in determining of the chiral coupling parameter L 3 , which at- 
tracts theoretical interest, extending beyond CHPT, for its direct relation to the gluon condensate 
and the constituent quark mass flUEl. The neutral K e4 branching ratio is mainly sensitive to L 3 
and very little to L 5 and L 9 lfl9l . Using this dependence one gets 



-4.1 ±0.2) x 10" 



(9) 



Also, in CHPT the form factor f p depends linearly on L 3 with directly computable numerical 
constants. In addition, the form factor A g depends linearly on L 3 with numerical constants de- 
pendent on the well known pion decay constant F^. Using either f p or A g , we get values for L 3 
consistent within errors with the value of eqn. © but with five times larger uncertainties. The 
value of L 3 determined in this work is more accurate than theoretical estimates from CHPT fits 
based on previously available data 0. Our result is also compatible with the result of ref. JH|. 
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